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Abstract. We address the issue of anomalous image flux ratios see n in the double-im age gravitational lens JVAS B0218+357. 
From the multi-frequency observations presented in a recent study iMittal et alllool and several previous observations made 
by other authors, the anomaly is well-established in that t he image flux-densi ty ratio (A/B) decreases from 3.9 to 2.0 over 
the observed frequency range from 15 GHz to 1.65 GHz. In lMittal et all feOOfili . the authors investigated whether an interplay 
between a frequency-dependent structure of the background radio-source and a gradient in the relative image-magnification can 
explain away the anomaly. Insufficient shifts in the image centroids with frequency led them to discard the above effect as the 
cause of the anomaly. 

In this paper, we first take this analysis further by evaluating the combined effect of the background source extension and 
magnification gradients in the lens plane in more detail. This is done by making a direct use of the observed VLBI flux- 
distributions for each image to estimate the image flux-density ratios at different frequencies from a lens-model. As a result 
of this investigation, this mechanism does not account for the anomaly. Following this, we analyze the effects of mechanisms 
which are non-gravitational in nature on the image flux ratios in B0218+357. These are free-free absorption and scattering, and 
are assumed to occur under the hypothesis of a molecular cloud residing in the lens galaxy along the line-of-sight to image 
A. We show that free-free absorption due to an Hn region covering the entire structure of image A at 1.65 GHz can explain 
the image flux ratio anomaly. We also discuss whether Hn regions with physical parameters as derived from our analysis are 
consistent with those observed in Galactic and extragalactic Hn regions. 



1. Introduction 

The double-image gravitational lens B0218+357 was discov- 
ered in the Jodrell- VLA 1 Astrq metric Survey (JVAS) of ra- 
dio sources JPatnaik et al.l 1 19921) . It has the smallest angular 
image-separation (~ 330 mas) amongst the known galactic- 
scale lens systems. The lensed object is a powerful radio galaxy 
(blazar) at a redshift of 0.944 JCohen et al.ll2003l) with a typi- 
cal core-jet morphology and a frequency-dependent structure. 
Further, it has a variable radio emission and the time-delay be- 
tween variations in the images has been accurately measured by 
iBiggs et al 1 lll999l) . (10.5 ± 0.4) days, wh ich is consistent with 
the value of (lO.l^j g) days meas ured bv ICohen et alJ d2000l) . 
The lens is at a redshift of 0.684 jBrowne etaTIIl993T see also 
O'Dea et al. 1992) and is believed to be a spiral galaxy based on 
the small angular image-separation and radio absorption lines 
that have been measured for this system. A robust confirmation 
is provided by the recent HST-ACS i mage of this syste m with 
a very high resolution and sensitivity jYork et al.l2005l) . which 
clearly shows the two point images and the underlying spiral 
structure of the lens galaxy. 

High-resolution maps of B0218+357 using long baseline 
interferometers, such as the VLBA and the VLBI networks, 
reveal A and B to consist of two distinct components with sim- 
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ilar separations (> 1 mas) but with different relative shapes 
and orientations. The double features, 1 and 2, shown in Fig.^ 
are identified with the core-jet morphology of a fiat-spectrum 
radio source, with component 1, based on its high turn over 
frequency, as the core or the jet-base, and component 2 as a 
jet-component. The phenotypical core-jet picture of the back- 
ground s ource is more pron ounced in the 8.4 GHz global VLBI 
maps bv lBiggs et al.l J2003I) . in which the low-brightness emis- 
sion constituting the jet is seen to extend out to ~ 15 mas to 
20 mas from the core in both the images. iMittal et"al1 d2006l) 
detected another component in the VLBI hybrid-map of image 
A at 1.65 GHz [see Figure (8a) of their paper], which is sep- 
arated by 12 mas from the superposition of components 1 and 
2. This feature is designated as component 3 and has not been 
found to have a counterpart in image B. 

All the observed characteristics of this lens system, such 
as the lens-geometry and the image positions, can be recon- 
structed well using a simple lens-model, except the image flux- 
density ratios. The anomaly in the image flux-density ratios 
in B02 18+357 wa s addressed and discussed extensively by 
IMittal et all J2006t hereafter M06). It does not violate the cusp 
or fold relations observable only in four-image systems, as is 
the case with the majority of other systems which show im- 
age flux ratio anomalies. Instead there is a steady decline in 
the image flux-density ratio (A/B) from 4 to 2 with decreas- 
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Fig. 1. VLBI 15-GHz map of B0218+357 for image A (right) 
and B (left), plotted with a restoring beam of 0.5 mas. 



ing radio frequencies from ~ 20 GHz to 1 GHz, which violates 
the achromatic nature of the phenomenon of GL. In M06, the 
authors took into account the combined effect of an extended 
background object with a frequency-dependent structure and a 
strong gradient in the relative image-magnification over the im- 
age plane, both of which exist in the B0218+357 lens system. 
The technique of inverse phase-referencing was used to estab- 
lish the positions of image-centroids at five radio frequencies. 
The results of their analysis led them to discard the gradient in 
the image-magnification ratio across the images as a cause of 
the flux ratio anomaly in B0218+357. 

In the present work, we extend this investigation but also 
seek other explanations, of entirely different origins, for the 
flux ratio anomaly. The intervening matter along the lines-of- 
sight to the gravitationally lensed images can, through non- 
gravitational (electromagnetic) effects, produce deviations in 
the image properties. Emission from the background source 
can be absorbed and/or scattered causing a change in the origi- 
nal radiation intensity. These mechanisms, in combination with 
the resolution of the observations, can affect the surface bright- 
ness of the lensed images differently and perturb the image 
flux-density ratio from its expected value. The most common 
physical processes that occur are extinction in the optical re- 
gion, and free-free absorption, scatter-broadening and Faraday 
rotation in the radio region. It is generally assumed that these 
mechanisms occur in the ISM of the lens galaxies and, when 
they occur, cause the lens to be no longer 'transparent'. It is 
to be noted that whichever mechanism is responsible for the 
flux ratio anomaly in B0218+357, it should explicitly pro- 
duce a frequency-dependent change in the image flux-density 
ratios, as observed. The physical quantities which appear in 
these above-mentioned astrophysical processes and which in- 
duce this v-dependence are the refractive index of the interven- 
ing plasma in the case of scatter-broadening and optical-depth 
in the case of free-free absorption. 

While the task of identifying the signatures of gravitational 
lensing becomes more difficult in the presence of above inter- 
fering effects, this can be exploited to our benefit by allow- 
ing us to probe the m edium of the intervening lens galaxy 
in detail (Bi ggs! 120051) . Based on radio propagation effects, 
VLBI techniques have been largely used to explore the lens 
galaxy in numerous lens systems. Such observations reveal 
a variety of constituents which otherwise are hard to detect, 



such as atomic and molecular species at differe nt redshifts 
dKanekar & Briggsl 120031 IChengalur etaD Il999l) and solar- 
mass objects that result in microlensi ng of the background radi- 
ation dKoopmans & de BruvnlEoOoT) . Radio polarization mea- 
surements help us to study the io nized gas fraction and large - 
scale magnetic fields in galaxies dNarasimha & Chitrell2004l) . 
These studies are interesting also for understanding the red- 
shift evolution of elemental abundances and large-scale mag- 
netic fields. 

There is a substantial body of evidence which indicates 
that the ISM in front of image A is rich in gas and dust. 
iFalco et all 1 1999b investigated 37 differential extinction curves 
in 23 gravitational lens galaxies and found that B0218+357 
and one other system (PKS 1830-211) in their sample have 
exceptionally high differential extinctions between the images 
of the lensed object (see also iMunoz et al.ll2004l) . According 
to their measurements, the selective extinction for image A is 
higher than for image B by AEb-v = 0.9 mag and the total 
(selective) extinction along the lines-of-sight to both the im- 
ages is equal to Eb-v = 0.6. This result is not surprising as 
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amounts of molecular gas and Hi in the lens galaxy. However, 
a strong relative extinction in image A additionally implies that 
the molecular cloud, which is associated with these molecu- 
lar absorption lines, lies in front of image A. This is in agree- 
ment with the observations of Hi absorpti on in B02 18+357 us - 
ing VLBI with a resolution of 80 mas bv lCarilli etaD J2000l) . 
They found that the dominant contribution to the Hi 2 1 cm line 
comes from the south- west component (image A ). Similar find- 
ings were achieved by iMenten & Reidl dl996l) who observed 
the formaldehyde (H2CO) absorption lines with the VLA at 
14.1 GHz and 8.6 GHz (the VLA provides an angular reso- 
lution of < 0.2 arcsec at these frequencies). They showed that 
the absorption arises solely due to image A and derived an up- 
per limit on the optical depth for image B, a factor three times 
smaller than that calculated for imag e A. This a l so len ds sup- 
port to the explanation proffered by lYork et all d2005l) in or- 
der to account for the mismatch in the image separation be- 
tween the radio and the optical measurements. Dust obscura- 
tion of image A can also explain why image B is observed to 
be brighter than image A at optical wavelengths (Lfackson et alJ 
l2000HLehar e7aTll2000l) while the opposite is true at all radio 
wavelengths. 

The goal of this paper is multifold. First, we carry our inves- 
tigation along the lines of argument similar to those presented 
in M06 further. The question that motivated additional investi- 
gation of the frequency-dependent structure of the B0218+357 
images is whether the magnification at the centroid position, 
li(x c ), gives a good estimate of the average magnification (to 
be defined in Section |5J suffered by an extended object. For 
a simple source structure this may reasonably and intuitively 
be assumed but for complicated source structures involving 
steep magnification gradients in the image plane, the assump- 
tion does not hold any longer. 
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Second, we investigate the effects of non-gravitational pro- 
cesses on the image flux-densities in B0218+357 in the hope 
of solving the long-standing problem of inconsistent image- 
magnification ratios. The main focus of the present work is to 
explore these non-gravitational effects. The layout of the paper 
is as follows. In Section[5]we give an introdu ction of the lens 
mode ls derived from the LensClean algorithm JWucknitz et alJ 
120041) . which we use throughout the paper. Here we also cal- 
culate the magnification-weighted image flux-density ratios in 
detail to assess the contribution of this effect to the flux-ratio 
anomaly. The propagation effects are described in Section 
which covers the free-free absorption mechanism (Section lXTt 
and refractive scattering (Section 13.21 . Finally, we present a 
discussion on the results obtained from the above analyses in 
Section @] and conclude our work on the flux ratio anomaly 
in B0218+357 in Section|5] All these investigations were car- 
ried out based on the multi-frequency VLBI observations of 
this lens system, which are presented in M06. Throughout 
the paper we adopt the flat ACDM concordance cosmological 
model for our calculations, with a Hubble parameter of Ho = 
71 km s _1 Mpc -1 and a cosmological constant of Qa = 0.73. 



2. Derivation of magnification-weighted image 
flux-density ratios using LensClean models 

The VLBI images (see M06) show that, whereas at 15.35 GHz 
the emission is dominated by the compact subcomponents 
with a separation of ~ 1 .4 mas, at lower frequencies (such as 
1.65 GHz and 2.25 GHz) a considerable amount of low bright- 
ness emission extends out to regions where the relative mag- 
nification is significantly different. Thus it may be insufficient 
to simply consider magnifications at the centroid positions as 
was the approach adopted in M06. The (true) resultant average 
magnification of an image is the integral of the (background 
source) magnification- weighted intensity over the image area, 



H v = 



J F v (x)d 2 x 
J [F v (x)/fi(x)] d 2 x ' 



(1) 



where F v is the measured flux component 2 at position x, fi(x) is 
the lens magnification from a model corresponding to this com- 
ponent and the integral is carried over the entire image area. 
The relative image-magnification is then simply 



V-R* = — 



2.1. Lens model 



(2) 



To use the above recipe for estimating the relative image- 
magnification demands the knowledge of the lens mass- 
distribution so that the magnification at all image points can 
be derived. The lens model used for the subsequent analysis 
is a singular elliptical potential with the mass-radius slope, B, 



2 The flux-density at a position x in the image plane, in principle, 
is not defined but in the following 'positions' will be identified with 
pixels that contain flux-densities. 



Table 1. The lens model parameters derived from the 
LensClean algorithm. The lens position coordinates are given 
relative to Al with right ascension increasing to the left (east) 
and declination increasing to the north. 



P 


xoi 
(mas) 


X02 
(mas) 


e\ 


ei 


ft 
(mas) 


1 

1.063 


255.214 
255.212 


117.193 
118.522 


0.0057 
0.0179 


-0.0494 
-0.0410 


163.269 
169.020 



fixed to 1 to obtain an isothermal profile (SIEP). The elliptical 
iso-potential form is given by 

2-jS 

if, (3) 



where 



(l +e) 2 (1 _ e) 2 



(4) 



Such a model is specified by five parameters, the lens position, 
xoi and X02, the ellipticity of the iso-potential contours of the 
lens, £\ and en, and the lens strength, ^o- The coordinate sys- 
tem adopted for the calculations is such in which x\ and xn are 
aligned with the major and minor axes of the ellipse [the lo- 
cus of points obtained with constant u in Eq. @] and the lens 
centre coincides with the origin. These parameters are derived 
from the LensClean algorithm and given in Tabled 

The isothermality can be disturbed by varying the value of 
B to values around 1 ; the best fitting non-isothermal value of B 
is given in Table ^ However, the procedure employed for this 
was not completely self-consistent. The lens position was con- 
strained to that determined from LensClean for an isothermal 
model and the other lens-parameters, including B + 1, were 
derived using the following VLBI constraints: the (Al - BY), 
(Al - A2) and the (B\ - B2) separations, and the image flux- 
density ratio at 15.35 GHz. 

There are two main motivations behind including non- 
isothermality that come from both multi-frequency observa- 
tions described in M06 as well as previous observations. Firstly, 
using the SIEP lens model on these data, there is an evident 
4 cr to 5 cr discrepancy between the observed and the modelled 
B2 — B\ component separation along right ascension. This re- 
duces to within 1 cr on applying the Singular Non-Isothermal 
Elliptical Potential (SNIEP). The second clue is obtai ned from 
a qualitative comparison made by Big gs et alJ 112003) between 
the CLEAN 8.4 GHz maps of images A and B back-projected 
into the source plane, wherein the B jet seems to be more elon- 
gated or stretched than the A jet by about 10 %. This dissimi- 
larity can be resolved only by invoking a different mass-radius 
profile, albeit the departure from isothermality needed to ac- 
count for this effect is very small. 

2.2. Method 

The model-predicted relative image-magnifications at different 
frequencies are calculated on the basis of the observed flux dis- 
tribution of either image A or image B (termed the primary 
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image). The flux distribution of the second image can be then 
derived using either of the lens models (SIEP or SNIEP). In 
this way two ratios can be derived by using information ob- 
tained from observations at each frequency. The first value is 
obtained by taking the ratio of the observed flux-density of A 
to the modelled flux-density of B and the second value is ob- 
tained by replacing the image A with image B and vice versa. 

The most tedious step in deriving the image magnifications 
is the inversion of the lens equation to solve for the image po- 
sitions for a given source position. To accomplish this, for ev- 
ery component in the primary image the corresponding com- 
ponents in the source plane were located. The secondary im- 
age components were determined by solving for the zeroes of 
the resulting equation. For the SIEP lens model, this can be 
achieved analytically but for the SNIEP model, since the result- 
ing equation is no longer a polynomial, the Newton-Raphson 
method was employed (for both purposes we used routines 
from the GNU Scientific Library, GSL). 

2.3. Results of detailed-averaging 

The results from the detailed analysis of the radio maps of the 
two images of B0218+357, A and B, and from applying the 
lens models are shown in Fig. |2] The line parallel to the x- 
axis is the image flux-density ratio of a background point-like 
source at the position of Al (or Bl) at 15.35 GHz, depicting 
the frequency-independent nature of the gravitational lensing 
effect. The model-predicted ratios for both isothermal and non- 
isothermal mass-radius profiles, obtained from using the ob- 
served brightness distributions of image A and image B sep- 
arately, are shown as square or triangle symbols. Also shown 
as 'stars' are the image flux-density ratios calculated from the 
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Fig. 2. Observed (star symbols) and model-predicted image 
flux-density ratios versus frequency. The open and filled sym- 
bols correspond to the isothermal and non-isothermal fitting, 
respectively. The square symbols are used when image A is 
taken as the primary image and the triangle symbols are used 
when image B is taken as the primary image. 



multi-frequency observations (M06), indicating clearly the fre- 
quency trend in the observed image flux ratios in B0218+357. 
As can be seen, for all four model-predictions the ratio remains 
the same to within < 4 % except for 1.65 GHz using A as 
the primary image, where the ratio differs from a point-source 
model by about 8 % for the isothermal model and 10 % for the 
non-isothermal model. Furthermore, the change in the ratio due 
to the shift in the centroid position in image A at 1.65 GHz, as 
derived from the phase-reference observations in M06, is less 
than a percent. 

Thus, our study of the effect of the interaction of frequency- 
dependent source structure with macro-model magnification 
gradient confirms that this cannot be the cause of the ob- 
served v-dependent image flux-density ratio in B0218+357. 
This analysis does show that the relative image-magnification 
estimated from applying the entire structure of image A at 
1 .65 GHz is significantly different from when the image- A cen- 
troid position at the same frequency is directly applied to the 
model. It is to be noted that the above 8 % to 10 % effect 
on the ratio is compatible with the values of relative image- 
magnification in the direction of the subcomponent newly- 
identified at 1.65 GHz only in image A, component 3. This is 
not surprising since the peak intensity of component 3 is com- 
parable to that of components 1 & 2, and has a non-negligible 
effect on the centroid position in image A at this frequency. 
However, the change in the ratio is in a direction in which 
the relative image-magnification increases, opposite to the de- 
clining trend observed at decreasing frequencies. Finally, this 
study confirms that there can be departures from the general no- 
tion that the phenomenon of gravitational lensing preserves the 
spectrum of the background source in all the images. These de- 
partures can be seen as small deviations of the relative image- 
magnifications at any given frequency relative to the neigh- 
bouring frequencies in Fig. [2] This has been shown based on the 
frequency-dependent radio structure of the background source 
in B0218+357 and without invoking any other external mech- 
anisms, some of which are described in Section|3] 

3. Propagation effects 

In this section the effects of free-free absorption and scattering 
on the properties of the B0218+357 images are investigated. 
Both these mechanisms require regions of plasma along the 
lines-of-sight to the image in consideration. The presence of 
a molecular cloud in front of image A provides an easy solu- 
tion to this requirement as molecular clouds harbour sites of 
recent star formation, which through photoionization build up 
regions of ionized hydrogen around them. 

3.1. Free-free absorption 

We begin with the assumption that the radiation from image B 
after lensing by the spiral galaxy is unperturbed by any other 
physical process. To investigate whether image A suffers from 
free-free absorption, we envisage a giant Hn region, embed- 
ded in the molecular cloud presented above, directly in front of 
the image. Hn regions are formed around hot O- or early B- 
type stars that emit uv-photons capable of ionizing the hydro- 
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gen clouds out of which these stars form in the first place. The 
delicate balance between the heating and cooling mechanisms 
within such a photo-ionized plasma, which comprises mainly 
hydrogen atoms, results in an equilibrium electron temperature 
of ~ 10 4 K. We further assume that this Hn region covers the 
entire image structure at the lowest frequency, 1 .65 GHz. This 
is the simplest scenario which can be considered in that it may 
well be modified to one consisting of several Hn regions re- 
siding within the same molecular cloud but at different spatial 
locations, which in projection cover the entire image. In the 
following we assume that the Hn region has a homogeneous 
distribution of electron density and a temperature in the range 
4000 K to 10 000 K 3 . 

For Hn plasmas a state of Local Thermal Equilibrium (LTE) 
is a good assumption. This is because the thermal equilibrium is 
well maintained in such regions as any perturbation that enters 
into the system is quickly redistributed through Coulomb colli- 
sions amongst all the elements of the plasma. This validates the 
Kirchhoff law, so that the free-free absor ption coefficient in the 
Rayleigh-Jeans reg ime is given by (see iRvbicki & Lightmai] 
ll979UBrownll987l) 

a s = 3.19 xlO- 7 f-^S-f f-?M (— I " Pc _1 (5) 
Ucm- 3 / \10 4 K/ \GHz/ V V ' 

where n e is the electron density, T e is the electron temperature 
and v is the observed frequency scaled by a factor of (1 + z\), 
where z\ is the redshift of the lens galaxy. In deriving Eq. l|5}, 
a pure hydrogen cloud has been assumed and thus, n e = ti\ and 
Z = 1, where «; is the ion density and Z is the ionic charge. 
The intensity of background radiation that passes through an 
Hn region is modified according to the law of radiative transfer 
to 

/ff(v,T) = / (v)^ T|r(v) , (6) 

where 1$ and / are the modified and background intensities, re- 
spectively, and Tff is the optical depth equal to the integral of the 
free-free absorption coefficient over the path length through the 
plasma. Under the assumption of homogeneity in electron den- 
sity and temperature, the optical depth is equal to a$L, where 
L is the depth of the cloud. In the above, the source function, 
which under the condition of LTE is equal to the Planck black- 
body spectrum, has been ignored since a source temperature of 
~ 10 4 K is below the minimum brightness temperature that can 
be detected using VLBI. Introducing the physical parameter, 
Emission Measure {EM = n^L), the optical depth is 

\10 4 K/ \cm- 6 pc/ I GHz/ V ' 

3.1 .1 . Parameter estimation of the Hu region 

In order to determine the parameters of the Hn region, we 
assume the spectrum of image B to reflect the true source 
spectrum. Then, using the isothermal lens model described in 

3 These values span the entire range of temperatures estimated 
for other galacti c and extragalactic Hn regions (see Chapter 3 of 
IOsterbroclJl974l) . 



Table 2. The parameters derived for the Hn region in front of 
image A. T is the temperature and EM is the emission measure 
of the Hn region. Given also are the various combinations of the 
electron density and the depth of the cloud for given emission 
measures. 



T EM 




L 


«e 


(K) (cm" 6 pc) 


(pc) 


(mas) 


(cm" 3 ) 


10 4 (1.8 ±0.3) XlO 7 


1 


0.15 


4243 




10 


1.5 


1342 




100 


14 


424 




200 


28 


300 


4 000 (5.3 ±0.9) xlO 6 


1 


0.15 


2302 




10 


1.5 


728 




100 


14 


230 




200 


28 


163 



Section 12.11 Q3 = 1 in Table Q and using the technique de- 
scribed in Section 12.21 the true spectrum of image A can be 
determined. Further, with the help of Eq. Q and from the 
knowledge of the true and the modified spectrum of image A, 
Eq. lHJ can be fitted to the observed flux-densities at different 
frequencies to determine the best fitting values of the plasma 
parameters. Unfortunately, the optical depth is degenerate to 
the following parameters which describe the Hn region: the 
electron density, the temperature and the depth of the cloud. 
But since the range of temperatures from observations of other 
Hn regions is known to be narrow, by fixing the temperature 
to two extreme values bracketing this range, the combination 
of the two remaining parameters, the emission measure, can be 
uniquely defined. The aim of this analysis is to find out whether 
the free-free absorption curve can be fitted at all to the observed 
spectrum of image A and also if the resulting values of EM are 
physically meaningful. 

The parameters of the hypothesized Hn region were esti- 
mated by, first, approximating the spectra of image B and the 
modelled image A by a synchrotron power-law, 

F A (v) oc v ~" ; F fl (v) oc v - b , (8) 

where Fa(v) is the image A flux-density and Fb(v) is the image 
B flux-density, and a = 0.153 ± 0.018 and b = 0.147 ± 0.022 
are the power law indices fitted to their spectra. We note that the 
spectral indices for image A and image B assume slightly dif- 
ferent values, in contradiction to the expectation of a constant 
magnification ratio for a point source from the model. This is 
due to the small differences in the modelled relative image- 
magnifications at varying frequencies, which in turn arise due 
to a frequency-dependent source structure (see Section 12.31 . 
Substituting Fa(v) for I (v) in Eq. (|6), the best-fitting value 
for EM can be calculated using the x 2 minimization method 
to minimize the difference between the observed and the mod- 
elled image A flux-densities. 

Shown in Fig. [3] are the flux densities of image A (open 
boxes) modelled using the observed flux densities of image 
B (crosses). The fitted power-law spectra are shown in dotted 
lines. The free-free absorption curve (solid curve) is fitted to 
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Table 3. Image flux-density ratios versus frequency as ob- 
tained from the observations (second column), the free-free ab- 
sorption model (third column) and the isothermal lens-model 
(fourth column). 



4 6 
Frequency (GHz) 

Fig. 3. The free-free absorption curve (solid curve) fitted to the 
observed spectrum of image A (filled boxes). The modelled 
flux-densities of image A (open boxes) follow the measured 
flux-densities of image B (crosses) and their spectra (dotted 
lines) have nearly the same power law index. 



the observed flux-densities of image A (filled boxes). The var- 
ious parameters, including the electron densities for different 
values of L, estimated for two values of temperatures are given 
in Table |2 The values of EM, although quite large, are con- 
sistent with those measured in giant Galactic Hn regions, and 
the estimations of « e and L are also in good agreement with 
those observed for Galactic and extragalactic Hn regions. In the 
present context the most meaningful combination of « e and L 
is the last entry in Table[2]for both T values. This is because we 
have assumed that the Hn region covers all of image A, which 
at 1.65 GHz (the lowest frequency) has a deconvolved size of 
~ 28 mas which translates into a physical size of about 200 pc 
at the redshift of the lens galaxy. 

Given in Table [3] is a comparison between the image flux- 
density ratios obtained from the observations (second col- 
umn) and those obtained after applying the free-free absorp- 
tion model (third column), using Eq. 0. There is good agree- 
ment between these values and it is seen that the absorption can 
explain the observed trend very well at all the frequencies ex- 
cept at 4.96 GHz where the two values differ slightly, as is also 
visible from Fig. [3] As a further comparison, the ratios calcu- 
lated from applying the SIEP model in Section l2~T1 for the case 
where image-B is the primary image, are also given (the fourth 
column). 

For the above analysis even though, as already mentioned, 
we assumed a complete and uniform coverage, we tried to put 
constraints on the position and depth of the Hn region relative 
to the brightness distribution of image A at different frequen- 
cies. This turned out to be difficult due to an insufficient num- 
ber of constraints available either from the observations used 
in this work (M06) or from the molecular (and atomic) line ob- 
servations of this system by other authors. Even so, the maps 
of image A from these observations are strongly suggestive of 
a part of the absorption screen being directly in front of the 
centre-most region containing components 1 and 2. We arrived 



v (GHz) 


pobs /pobs 
r A ' r B 


pit /pobs 


pmod /pobs 
r A /r B 


1.65 


2.07±0.06 


1.99 


3.90 


2.25 


2.71+0.11 


2.73 


3.80 


4.96 


3.01+0.06 


3.46 


3.76 


8.40 


3.41+0.15 


3.64 


3.82 


15.35 


4.01+0.15 


4.01 


3.84 



at this conclusion by mapping the inner subregions of image 
B (centered on Bl), 2 mas by 1 mas in size, into equivalent 
regions in image A at all five frequencies. A comparison be- 
tween the modelled spectrum of image A and the observed one 
showed that while there is a good match between the modelled 
and the observed flux-densities at 15.35 GHz, there is a clear 
indication of strong absorption at lower frequencies and that the 
fraction of the absorbed continuum increases with decreasing 
frequency. 

From the above discussion, it can be concluded that the 
free-free absorption hypothesis is capable of reproducing the 
observed spectrum of image A and, thereby, of solving the im- 
age flux-density ratio anomaly in B0218+357. Furthermore, 
the values of the emission measure resulting from the fit for 
two extreme electron temperatures are quite reasonable in that 
similar values have been measured for Galactic and extragalac- 
tic Hn regions, lending further support to the hypothesis. 

3.2. Scattering 

Extragalactic sources are often found to be scatter-broadened 
by an intervening screen of turbulent plasma. The fluctuations 
in the plasma density induce variations in the refractive index, 
and in turn scatter the background radiation. This gives rise 
to two effects, scintillation and image- or scatter-broadening. 
There have been claims of scatter-broadening se en in grav- 
itational lenses as well, s uch as PKS 1 830-21 1 jjones et alJ 
Il996l iGuirado et alJll999|) B 193 3+503 jMarlow et al.ll 19991) 
B02 1 8+357 jBiggsetall l2003t). B0128+43 7 jBiggs et alJ 
120041) and PMN ,T1838-3427 (fwinn et 3112001 . 



Scattering by an infinite screen with homogeneous fluctua- 
tions in the electron density does not result in any changes in 
the flux-density of the background source. Based on symmetry 
arguments it can be proven that, in such a case, as much flux is 
scattered away as toward the line-of-sight to the observer, re- 
sulting in the conservation of flux-density. But this is no longer 
the case if there are variations in the statistical properties of the 
scattering screen. This includes the case where the scattering 
material appears in front of the scattered image only in parts, 
separated by 'holes' comprising neutral matter. In particular, 
if the scattering screen is truncated over the size of the scat- 
tered background i mage, it will lead to a n attenuation of the 
source flux-density JCordes & Eaziol200ll) . Thus, variations in 
the lensed image flux-densities because of an intervening scat- 



Rupal Mittal et al.: Free-free absorption in the gravitational lens JVAS B0218+357 



7 



tering medium can occur only when there are discontinuities or 
variations in the scattering strength. 

Scattering studies and measurements have been popular for 
over two decades now, and the basic physics behin d the scat- 
tering mechanism is quite well -established (e.g. iRickettll 1 99(1 
iNaravarl 1 99llGoodmanl 19971) . The regime of relevance in the 
present context is refractive scattering, as opposed to diffrac- 
tive scattering. Refractive scattering is due to large-scale in- 
homogeneities in the medium with length scales, r le f, known 
as the refractive scale. This scale is much greater than what is 
known as the diffractive length-scale, r^ff, which represents the 
transverse spatial-scale for which the root-mean-square phase 
change due to the intervening plasma blobs is equal to 1 rad. 
For refractive scattering, 6 r = r K i/D\ corresponds to the size of 
the scatter-broadened image of a point source projected back 
on the scattering screen (which is assumed to be in the lens 
plane), where D\ is the angular diameter distance from the 
observer to the scattering screen. For a source with a non- 
negligible size in comparison with the angular broadening size, 
the Gaussian-equivalent angular size of the image as measured 
is the quadrature sum of the intrinsic source-size and the re- 
fractive length scale. The latter scales as v~ 2 - 2 [see Eq. J12H . 
which introduces the v-dependence of the scattering observ- 
ables. Diffractive scattering, on the other hand requires suffi- 
ciently compact sources, 6 < r^/Di, which for extragalactic 
scattering screens is rarely ever fulfilled, and therefore is not 
considered here. 

To probe the relevance of scattering, the first issue that is 
addressed in the following is whether image A is affected by 
scattering at all. From an observational standpoint, the quanti- 
ties which provide evidence for scattering are the image-sizes, 
in particular their dependencies on frequency, and temporal 
variations in the image flux-density. The latter can be ruled out 
since the image flux-density ratios measured for B02 18+357 
are independent of the time of observations, apart from the 
5 % to 10 % level fluctuations arising due to the source vari- 
ability and the time-delay between the images. To determine 
whether image A is scatter-broadened relative to image B (or 
vice-versa), we fitted two-dimensional elliptical Gaussians to 
the smallest detectable components in both the images at all 
frequencies. Shown in Fig. @] are the equivalent circular radii, 
Ra and Rb, of the ellipses fitted to these components plotted as 
a function of frequency. On fitting power-law curves through 
these data points ~ y~*° J ', where the power law indices, k a and 
kb, correspond to the curves through A and B respectively, it is 
found that 

k a = 1.40 ±0.11 ; k b = l.\l ±0.05. (9) 

It should be pointed out that the scattered images produced 
by a homogeneous scattering screen, in general, have circular 
shapes due to isotropic scattering. Image A, on the contrary 
is elongated in roughly the same direction as predicted by the 
magnification matrix. In other words, the tangential eigen di- 
rection due to lensing is preserved. This is not an argument 
against scattering, however, because the phenomenon of grav- 
itational lensing affects the scattering angle in the same way 
as the background source (see Appendix lAl. i.e. the scattering 
angle also gets magnified in the tangential direction. 



The presence of two images of the same background 
source, if only one is assumed to be affected by scattering, can 
be used to advantage through determination of the true source 
size using the unaffected image. To determine the effects of 
scattering on the B0218+357 image properties, similar argu- 
ments are followed as those presented in Section |3~T1 in that 
it is assumed that image B is unaffected by scattering or any 
other mechanism such as free-free absorption or milli-lensing 
(see Section^). Then, the size of a component measured at any 
frequency in image B projected back to the source plane us- 
ing a model will correspond to the intrinsic size of this compo- 
nent. From this, its equivalent size in image A can be calculated 
and compared to the observed size to yield a value for the re- 
fractive scale at that frequency. But this one-to-one correspon- 
dence can, from these observations, only be made at 15 GHz 
since at lower frequencies neither of the components 1 and 2 
can be resolved, nor has any other distinct component been un- 
ambiguously detected in both the images. The third compo- 
nent in image A does not add any new information since it can 
either be that it is intrinsic to the source structure but due to 
insufficient resolution remains undetected in image B, or that 
it is produced by some other astrophysical mechanism. Thus, 
translating image B component-sizes to the true source sizes 
is not a reliable approach for analyzing the scattering hypoth- 
esis. Therefore, even though the similarity between the fitted 
power-law indices for both the images is highly suggestive of 
no, or hardly any, scattering in front of image A, the possibil- 
ity cannot be ruled out completely. Accordingly, the scattering 
supposition in front of image A is pursued a little further in the 
hope to derive some meaningful limits. 

Before deriving limits to # r , it should be recalled that the 
observed component-sizes are magnified by the lens galaxy. 
For a rough translation of a measured size in the image plane 
to its true size in the source plane, the following recipe is used. 
Consider a small elliptical background object with a' as the ma- 
jor axis and V as the minor axis. Let the major axis be aligned 
with the direction of the tangential magnification. Then, for 
an isothermal mass profile, the source will be lensed into an 
elliptical image of size J ijia') x b', where is the tangen- 
tial magnification (the radial magnification for an isothermal 
mass profile is unity). Equivalently, if an image-component 
has a size of yaxb, its equivalent size in the source plane 
is J(a X b)/fi. For the current analysis, we use approximate 
values of 2 and 0.5 for the tangential magnification and de- 
magnification for images A and B, respectively. Consequently, 
the estimated sizes of the source-component 1, from the equiv- 
alent circular sizes of Al and Bl at 15 GHz, are: 

SA = 0.25 ± 0.03 mas ; 9 SB = 0.28 ± 0.03 mas . (10) 

Thus, there is a good agreement between the observations and 
the model-predictions at 15 GHz. But this is not surprising, 
since it is generally accepted that 15 GHz observations of this 
lens system are devoid of any discrepancies. The lens models 
discussed in Section |2~T1 are derived assuming that the image 
flux-densities and positions observed at 15 GHz are free from 
any intervening astrophysical mechanisms since most interstel- 
lar effects go as v~ 2 . 
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Fig. 4. Log-log plots of the geometric means of the axes of the elliptical Gaussians fitted to the smallest components in image 
A (right panel) and image B (left panel) versus frequency. Note that the scales on the y-axes are different for the two images. 

Table 4. The scattering angular-sizes as a function of frequency for the Scattering Measure, S M = 2.12 kpc irT 20 / 3 (T) and 
SM = 0.76 kpc irT 20 / 3 ('II'), the quadrature difference of t9 SA and 6^ U y the component sizes in the image plane, 9a and 9b, and 
the source plane, 9 $A and t9 SB [see Eq. i ll It and Eq. J12H 



v (GHz) 


9 r (mas) 


9 r (mas) 


6>a (mas) 


6» Sa (mas) 




8b (mas) 


6 SB (mas) 




I 


II 
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II 












1.65 


32.40 


17.47 


7.29 


3.93 


10.31+0.01 


7.29 


6.14 


4.34+0.37 


6.14 


2.25 


16.38 


8.83 


3.70 


1.99 


8.20+0.08 


5.80 


5.45 


2.60±0.01 


3.68 


4.96 


2.88 


1.55 


0.65 


0.35 


2.81+0.05 


1.99 


3.84 


1.05+0.01 


1.48 


8.40 


0.90 


0.49 


0.20 


0.11 


1.00+0.03 


0.71 


0.70 


0.40+0.03 


0.57 


15.35 


0.24 


0.13 


0.05 


0.03 


0.36+0.01 


0.25 


0.25 


0.20+0.04 


0.28 




On the other hand, assuming that the component size of im- 
age A fitted at 1.65 GHz is solely due to scattering, an upper 
limit of <9 r - 9 Sa = 7.29 mas is derived. It is to be noted that 
this is the observed scattering angle and depends upon the lo- 
cation of the scattering screen with respect to the source and 
the observer. The true scattering angle is 



A 
A 



r 



where A 1S the angular diameter distance between the scat- 
tering screen (the lens galaxy) and the source. The quantity 
which is sought is the Scattering Measure (SM, in units of 
kpc rrT 20 / 3 ), which is the integral of the strength of the tur- 
bulence in the electron density distribution along the line-of- 
sight. Following the mathematical prescription given bv lWalkerl 
Jl998lEoOll) for a homogeneous scattering screen, the relation 
between the refractive scale for a point source ( < 9 r ) and S M 
is 



SM = 



\ mas / I GHz/ 



5/3 



kpc m 



■20/3 



(12) 



where v is the radiation frequency at the redshift of the lens 
galaxy. Thus as pointed out earlier in this section, the scatter- 
ing angle corresponding to the refractive scale, 9 r [which is re- 
lated to 9 r through Eq. JTTi l. scales as v~ 2,2 . Substituting 9 r and 



v = 1.65 GHz in Eq. O, a value of SM = 14.2 kpc rrT 20/3 is 
obtained. It should be pointed out that the above limit on 9 r is 
a secure upper limit since the intrinsic size of this component 
at 1.65 GHz is > 0. In fact, 9 Sa = 7.29 + 0.01 mas is close to 
t9 SB = 6.14 ± 0.52 mas. Therefore, the true intrinsic source-size 
corresponding to the smallest component measured in image A 
at 1.65 GHz must be close to t9 SA . This will decrease the scat- 



(11) tering measure. 



A second value of SM can be estimated by assuming 
a true correspondence between the observed components at 
1.65 GHz. Then, the quadratic difference in the back-projected 
sizes of the components can be attributed to the intervening 

scattering disk. In this case, 9 r = -^t9 2 A - t9 2 B = 3.93 mas. This 

leads to a value of SM = 5.1 kpc rrT 20 / 3 . 

Substituting these values of scattering measure into 
Eq. (I12> . values of c9 r can be determined at other frequencies 
as well and are given in Table |4] T corresponds to when 
the A component at 1.65 GHz is used and 'II' corresponds 
to when the back-projected component in B at 1.65 GHz is 
assumed to represent the intrinsic source-size. Also given in 
Table |4] are the sizes of the image components in the image 
plane, 9& and 9b, and projected back to the source plane, 6> Sa 
and # SB , and the quadrature difference of 9 SA and t9 2 n) . For the 
first method (T), 6 t and i9 SA can be directly compared and it is 
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clearly seen that except at 1 .65 GHz, there are strong inconsis- 
tencies at other frequencies. For the second method ('II'), the 
difference between 9 SA and 6*,( U) , in quadrature, should be com- 
pared with SB . Once again, except at 1.65 GHz the values don't 
match. It should be mentioned, nevertheless, that the two val- 
ues for S M derived in this section are comparable to the values 
of S M calculated f or other galactic and extragalactic scatter- 
broad ened sources JCordes & Lazioll2002l ISoangler & Cordesl 
fl99aiSnan g leretalJl986l). 

iBiggs et all J2003I) have suggested a much higher value of 
SM = 150 kpc nr 20/3 for B0218+357 which is atypical com- 
pared to those measured along typical lines of sight through the 
Galaxy but not inconsistent with the measurements in the direc- 
tion of the Galactic center. However, their value is clearly not 
in agreement with the upper limit derived from these observa- 
tions, given in Table |4] On the other hand, on account of free- 
free absorption, of which there is good evidence, it may well 
be that the estimated sizes in component A at lower frequen- 
cies are somewhat different from those due to scattering alone, 
which weakens any argument against the scattering hypothesis. 
In addition, it can also be that the scattering disk is such that 
it covers components 1 and 2 only, so that it has an observable 
effect on the size-measurements at 8.4 GHz. From comparing 
the back-projected sizes of components in A and B at 8.4 GHz, 
we find that a slightly bigger size of the image-A component 
requires a scattering measure of « 48 kpc trr 20 ' 3 . In such a 
case, the effect of scatter-broadening would neither contradict 
the size-measurements at 15 GHz, on account of the 9 r ~ v~ 2 - 2 
frequency scaling, nor the low-frequency size-measurements, 
on account of a large intrinsic source-size. Such a scenario will 
be compatible with both, the claim made by Bi ggs et alJ (12003) 
that image A is scatter-broadened at 8.4 GHz and also our mea- 
surements given in Table[4] 

To summarize, it seems that scattering is not a satisfactory 
explanation for the flux ratio anomaly seen in B02 18+357, even 
though it might contribute to scatter-broadening of the core- 
jet components in image A at higher frequencies. The reasons 
are two-fold. First, there is no strong evidence of v-dependent 
image-broadening for image A. This follows from the similar- 
ity of the v-dependencies of the component-sizes in both the 
images and under the assumption that image B does not suf- 
fer from any scattering. Moreover, the v-dependencies do not 
resemble the standard v~ 2 scattering law. Instead, the index- 
values (slightly bigger than one) indicate that the increase in 
image-sizes with wavelength is intrinsic, common to most syn- 
chrotron self-absorbed sources. Second, even though the scat- 
tering measures derived in the above analysis are within the 
observed range of scattering measures in other systems, the 
back-projected component-size in image A and the size of the 
scattering disk are not compatible with each other at all five 
frequencies. However, these conclusions may be drawn only if 
the scattering strength is assumed to be homogeneous through- 
out the extent of the image A. If the scattering screen is clumpy 
and anisotropic, the analysis is no longer as straight-forward as 
that presented here. 



4. Discussion 

Two propagation effects, free-free absorption and scattering, 
have been investigated in the context of the flux ratio anomaly 
in the lens system B0218+357. Both these ISM effects have an 
inverse-frequency-squared dependence. In free-free absorption 
it is the optical depth that scales as v -2 ' 1 , while in scattering it is 
the effective size of the scatter-broadened image that scales as 
y-2-2 w e assume that they occur in the interstellar-medium of 
the lens galaxy. In principle, these mechanisms can occur any- 
where along the line-of-sight but the presence of a molecular 
cloud at the redshift of the lens galaxy lends the first possibility 
a higher probability. The image B properties, the flux-densities 
and sizes at varying frequencies, are assumed to remain unaf- 
fected by these processes. 

To test the free-free absorption hypothesis, an Hn plasma in 
front of image A is envisaged which comprises free electrons 
that through Coulomb collisions absorb a part of the incoming 
image-A radiation, 7 . The image A flux-density is attenuated 
by a factor, e~ T{y \ which depends upon the optical depth (t) of 
the cloud, and the spectrum of image A is altered. Using the 
isothermal lens-model described in Section |2~T1 the observed 
flux-densities of image B are used to calculate the correspond- 
ing (unabsorbed) flux-densities of image A at all the frequen- 
cies. These values are used to fit the free-free absorption curve 
to the observed flux-densities of image A (Hig. |3j and the EM 
of the Hn region is constrained to within 13 % accuracy. This 
has been done for two values of the electron temperature, which 
is known to have a very narrow range from copious observa- 
tional data of other Hn regions. The resulting values of EM 
appear to be in the typical range observed in other galactic and 
extragalactic Hn regions. Thus, free-free absorption is an excel- 
lent candidate to explain the frequency-dependent image flux- 
density ratios. 

From the lens models described in Section |2~T1 image A 
lies at a projected distance of 2 kpc from the lens centre. Its 
line-of-sight may well pass through one of the spiral arms of 
the galaxy, which are known to harbour extensive star forming 
regions and are rich in population I objects, such as young blue 
stars surrounded by Hn regions. In fact observations of other 
late-type spiral galaxies, such as M 51 jFoet alJll98ll . have 
shown that molecular gas, in the form of huge complexes of gi- 
ant molecular clouds, is strongly confined to the spiral arms of a 
galaxy. Thus, it is possible that the image-A line-of-sight passes 
through such a region rich in gas and ionized plasmas (Hn 
regions). From the measurements of NH3 absorption lines, 
iHenkel et al.l J2005I) visualize the molecular absorber in front 
of image A to be elongated along a path with roughly constant 
galactocentric radius to reconcile with high source covering 
factors observed at mm- wavelengths. Elongated filament-like 
molecular clouds have been observed in the Milky Way, such 
as the Orion Giant Molecular Cloud (GMC) and the Monoceros 
GMC. For the free-free absorption hypothesis, it has been as- 
sumed that the Hn cloud covers image A completely at all the 
observed frequencies. At 1 .65 GHz, the total extent of image 
A in the tangential direction is ~ 28 mas which corresponds 
to a linear size of ~ 200 pc in the lens plane. Such giant (and 
supergiant) Hn regions, although not ubiquitous, have indeed 
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been observed, both as galactic and extragalactic. For example, 
the giant Hn region complex W49 in the Milky Way Galaxy 
has L = 150 pc and n e = 100 cm 3 , NGC 604 in M33 has 
L = 400 pc and n e < 60 cm 3 and NGC 547 1 in M101 has 
L = 800 pc and n e = 200 cm 3 dShieldslll990h . Thus, the sce- 
nario in which the flux-density of image A is perturbed by a 
giant Hii region embedded in a turbulent molecular complex 
via free-free absorption can be easily envisioned. In reality, it 
may well be that there are several Hn regions embedded within 
the same molecular cloud but at different spatial locations and 
in projection cover the entire image A. Lastly, the question 
arises whether such a GMC in front of image A can, by virtue 
of its large mass, result in the magnification of image A, and 
possibly also a change in its posit ion, differing from that pre- 
dicted by the current lens-model. ISolomon & SandeTsl Jl980l) 
estimated the mass-spectrum of the observed clouds in the 
Milky Way Galaxy and concluded that about 90 % of the mass 
of the molecular ISM is contained in GMCs with sizes larger 
than 20 pc and typical masses larger than 10 5 Me- On the other 
hand GMCs are not very centrally concentrated. Even though 
they might contain regions of high surface-density, such as 
compact vigorously star forming clumps, such scattered mass- 
condensations do not act as efficient lenses. Therefore we can 
ignore the gravitational effect of a plausible giant molecular 
complex on the properties of image A. 

As a separate issue, we investigated whether there is any ev- 
idence of contribution from scattering in the lens galaxy to the 
disagreement between the predicted and the observed image 
flux-density ratios. The scatteri ng hypothesi s has b een invoked 
for this lens system before by IBi ggs et all who com- 

pared the 8.4 GHz maps of the images back-projected on the 
source plane and concluded that the subcomponents Al and A2 
have bigger sizes than Bl and B2. This, they claimed, was due 
to scatter-broadening in image A due to a turbulent medium 
in the lens galaxy, which is known to be ion and gas rich. 
In Section l3~2l similar comparisons between the image sizes 
have been tried with the present multi-frequency observations 
by applying a simple but effective method of determining the 
equivalent circular sizes of the components in the source plane 
from their fitted values in the image plane (the last two columns 
given in Tabled. But a direct comparison of the source-sizes 
derived from image A with those from image B presents a dif- 
ficulty due to the ambiguity in the identification of the corre- 
sponding components in the images. Except at 15 GHz, due 
to different image magnifications it is not clear if the elliptical 
Gaussian fitted to an unresolved component in one image cor- 
responds to the superposition of the same set of components in 
the source as fitted in the other image. The 8.4 GHz maps that 
IBiggs et all J2003h used for their analysis were obtained from 
global VLBI observations of B0218+357 with an rms noise 
of 30 /uiy beam -1 . With this sensitivity, components 1 and 2 
can easily be resolved in both the images. The rms noise in 
the 8.4 GHz maps of B0218+357 obtained from observations 
used in this work, in comparison, is at least a factor 10 higher 
and the components are unresolvable in image A and just re- 
solvable in image B (this is because the component-separation 
along R.A. is larger in image B than in image A). The sizes 
of the components fitted to the images indicate, individually, a 



similar frequency-dependence, y~*« b , with k\, = 1.2 for image 
B, and a slightly stronger dependence with k. d — 1 .4 for image 
A. It is to be noted that while the scattering angle scales with 
the inverse of frequency-squared, the intrinsic size has a much 
less pronounced dependence on the frequency. Based on exten- 
sive observational evidence it is seen that the angular size of a 
synchro tron self-absorbe d source scales roughly as inverse fre- 
quency jMarscheJ 1 9771) . Thus, the increase in the component- 
sizes in both the images can be an intrinsic effect. On the other 
hand, the scattering medium considered in this analysis is as- 
sumed to have statistically homogenous properties. It may well 
be that this assumption is not valid and that either the scattering 
measure has spatial variations transverse to the line-of-sight or 
the scattering screen is patchy with clumps of ionized mate- 
rial distributed over the image plane. Depending upon the dis- 
tribution of ionized clumps, the image flux-density can either 
amplify or attenuate. In addition to this, effects from free-free 
absorption can further complicate the investigation of the scat- 
tering hypothesis by disturbing the flux-density distribution in 
image A, and thereby, altering the size-measurements derived 
from the elliptical Gaussian fits to the components. Therefore, 
even though there is no hard evidence for scattering found in 
our data, which dominates either the image sizes or the image 
flux-densities, we can not completely eliminate the possibility 
of a partially-covering scattering-screen. 

One other plausible mechanism, which can lead the im- 
age magnifications to deviate from values dictated by sim- 
ple macro-lensmodels, and which we have not addressed so 
far is milli-lensing or substructure due to the subhalos in the 
mass range from 10 6 M to 10 9 M Q . Such subhalos are pre- 
dicted by semi-analytic and numerical si mulations of galaxy 
formation based o n a CDM universe (e.g. iKlvpin etalj 1999; 
iMoore et al Jl999i) . The number of such substructures are over- 
predicted by almost an order of magnitude compared to what 
is observed around the Milky Way Galaxy. The anomalous 
image flux ratios in galaxy scale lens systems are thought to 
be the first direct evidence of such small-scale power (e.g. 
iMao & SchneideJl998tl5alal & Kochanekl2002l) . But for two- 
image lens systems substructure effects cannot be examined 
with ease because of a severe shortage of constraints. In prin- 
ciple, a distribution of CDM-subclumps around one or both 
the images with a certain range of masses can reproduce any 
observations comprising image-positions and flux-ratios. But 
such a solution is not unique because in such a case the num- 
ber of model parameters exceeds the number of constraints 
as obtained from observations. However, it is also predicted 
that the CDM mass-clumps around the lensed images pro- 
duce astrom etric shifts in the centroid of the image brightness- 
distribution JPobler & Keetonll2005l) . These shifts depend on 
the source structure (thus frequency) and can be easily detected 
with high-resolution VLBI observations. The peak-to-peak im- 
age separation in B0218+357 does show an anomalous shift- 
ing at lower frequencies by about ~ 2 mas, therefore, entail- 
ing a need to follow-up the substructure hypothesis. For future 
studies, one way of analyzing the substructure hypothesis in 
B0218+357 further is by fixing the macro lens-model to the one 
obtained from the LensClean algorithm and generating random 
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realizations of substructures JPalal & Kochanekll2002l) based 
on analytic approximations of substructure statistics. 
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Fig. 5. Slices made through the core-components (1 & 2) and 
component 3 in image A at 1.65 GHz at a resolution of 7 mas 
and 2.25 GHz at a resolution of 5 mas. 



The anomalous frequency-dependent peak-to-peak image 
separation might also find an explanation within the frame- 
work of free-free absorption. It may well be due to inhomo- 
geneities in the absorbing-screen, which results in the centroid 
of the brightness distribution in image A moving toward the 
south-east direction, thereby, reproducing the separation trend. 
An inhomogeneous absorbing-screen may also provide an ex- 
planation for the occurrence of component 3 in image A at 
1 .65 GHz (as an alternative to the explanation offered in M06, 
namely that component 3 is a distinct real part of the back- 
ground structure and does not appear in image B due to in- 
sufficient resolution). The VLBI maps of image A obtained 
from our data indicate a gradual appearance of the secondary 
maximum corresponding to component 3 with frequency. Even 
though it is visible most prominently at 1.65 GHz, there are 
hints that it exists even at 2.25 GHz. This is seen as an un- 
resolved shoulder in the intensity cut at 2.25 GHz in image A, 
shown in Fig. [3] along the direction containing the core-jet sub- 
components and component 3. It seems that the absolute inten- 
sity of this component at both frequencies is nearly the same 
but that the fractional intensity relative to the peak intensity 
(centered around the core-jet subcomponents) decreases from 
1.65 GHz to 2.25 GHz. This might be due to a lower value of 
the absorption coefficient in the region of component 3 rela- 
tive to the rest of the image area, which would have the effect 
of making the corresponding feature of the background source 
more prominent with decreasing frequencies. 

5. Conclusions 

We have shown that the image flux ratio anomaly in 
B0218+357 can be explained by invoking propagation effects, 
namely those arising due to free-free absorption, which is 
thought to occur in an ionized medium of the lens galaxy. It 



is believed that this is the first time that there has been of- 
fered a plausible explanation for the image flux ratio anomaly 
seen in B0218+357. This result has been achieved solely on 
the grounds of multi-frequency radio continuum observations 
presented in this work. The hypothesis receives support from 
radio and optical observations of molecular and atomic absorp- 
tion lines which indicate the presence of an exotic interstellar 
medium in the lens galaxy, especially in front of image A. 

A multi-frequency VLBI analysis of B0218+357 has led us 
to conclude that in order to investigate the causes of flux ratio 
anomalies evident in numerous lens systems, a multi-frequency 
approach is most yielding. Many lens syst ems are thought t o 
be affected by dust-extinction in the optical jFalco etalll999h . 
and by scintillation and scatter-broadening due to small-scale 
inhomogeneities in the ionize d component of th e ISM in the 
radio e.g. PMN J l 838-3427 JWinn etal.ll2004. B0128+437 
llBiggs etaj]l2004 and PKS 1830-21 1 Jduirado et aljfl999l 
I.Tones et alJll996h ~ The process of free-free absorption affects 
the spectra of the images causing them to be different from 
each other at low frequency, especially if the line-of-sight to an 
image passes through the galaxy centre or its spiral arms where 
the el ectron column den sities are large, such as in PMN J 1632- 
0033 JWinn et al.ll2003i) or, as has been shown in this work, 
JVAS B0218+357. Even milli-lensing, though gravitational in 
origin, can lead to v-dependent deviations in the properties of 
the lensed images. Previous CDM-substructure investigations 
have been made at a single frequency implying a frequency- 
independent effect of substructure on the image flux densi- 
ties and image positions. However, the magnitude of the sub- 
structure effect on image magnifications changes with varying 
source size (especially apparent in the radio) and position rel- 
ative to the Einstein radius of the perturber. Observations at 
several frequency bands at radio and optical wavelengths are 
essential to probe different processes, with different frequency- 
dependencies, acting on the image magnifications and to assess 
their contributions. 

Further, we note that in order to use image flux-density ra- 
tios measured at radio wavelengths as constraints for lens mod- 
elling, the most trust-worthy values are those measured at high 
frequencies (such as 15 GHz and above) since the propaga- 
tion effects at these frequencies are minimal. Also, any effects 
of substructure , such as different morphologies of jet images 
( lMetcaldl200l . become more easily visible at such frequen- 
cies, which serve as further useful constraints for calculating 
substructure masses. 
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Appendix A: Magnification of the scattering angle 

In the presence of gravitational lensing, caused by an inter- 
vening massive object between the source and the observer, 
the scattering mechanism does not occur independently. We 
present a simple argument to show that the scattering angle is 
affected by lensing in the same way as the shape and size of the 
background source. 

Point Source 




Fig. A.l. Scattering of a point source in the lens plane. 

Consider a one-dimensional case of a point source which is 
scattered by the medium of an intervening galaxy. Let the scat- 
tering size be 9. As the next step, we introduce the deflections 
contributed by gravitational lensing which results in (0/2 — a\) 
and (8/2 + a 2 ) as the total deflection at the edges of the lensed 
image, where a\ and a 2 are the angular deflections attributed to 
lensing. Then, the condition for the light rays from the lensed 
object to reach the observer is, 



da 

9obs = + (a 2 -ai) = 0+— fobs 

d0 



(A.l) 



where o b s is the observed angular extent of the lensed object. 
Hence 



= n0 



1 



V- = 



1 - da/dd ' 



(A.2) 



where \i is t he radial magnification produced by an axially sym- 
metric lens JSchneideretai]ll992l) . Thus, the observed image 
size is the scattering size scaled by the lens-magnification fac- 
tor. The result just derived does not apply only for radial mag- 
nification but is true in general. 
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